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Abstract 

Some proposals to account for the highest energy cosmic rays predict 

that they should point to their sources. We study the five highest energy 

events {E > 10^^ eV) and find they are all aligned with compact, radio-loud 

quasars. The probability that these alignments are coincidental is 0.005, given 

the accuracy of the position measurements and the rarity of such sources. The 

source quasars have redshifts between 0.3 and 2.2. If the correlation pointed 

out here is confirmed by further data, the primary must be a new hadron or 

one produced by a novel mechanism. 
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The nature and origin of the highest energy cosmic rays {E > 10^° eV) is one of the major 
questions in physics and astronomy. Energies up to 3.2 x 10^° eV corresponding to center- 
of-mass energies up to ^/s ~ 800 TeV, have been observed for the primary interaction with an 
atmospheric nucleon. The showers produced by these cosmic rays indicate that the primary 
is a hadron such as a proton or hght nucleus although a photon is not completely 

excluded. Astrophysical mechanisms to accelerate protons to energies of up to io^^~^^ eV 
have been identified [Q, but they require exceptional sites. In his pioneering analysis, Hillas 
P observed that the source could be a radio galaxy or quasar, and not much else, based on 
general considerations. 

The conundrum is that nucleons, nuclei, and photons of energy greater than about 5x 10^^ 
eV have a non-negligible scattering cross section from the cosmic background radiation 
(CBR), causing their energy to be reduced to this level if they travel far enough through 
the CBR. This is known as the Greisen, Zatsepin, Kuzmin (GZK) limit If a proton 
or photon arrives at Earth with an energy greater than 10^° eV, it is exceedingly unlikely 
to have originated further than 50 Mpc whereas suitable astrophysical acceleration 

sites are located at greater distances 0, and indeed none are found within the expected 
scattering cone of the higest energy event at less than the GZK distance [I7[. A heavy 



nucleus might survive the requisite distance and would have a larger scattering cone, 
but it is difficult to reconcile the observed flux with the < 10~^ Fe to p ratio corresponding 
to cosmic abundances. 

Proposals to resolve the puzzle range from positing superheavy relics - topological defects 
or heavy particles - whose decay produces nucleons and photons within the GZK distance 
|ll|], to positing new particles or mechanisms which evade the GZK bound [p|, p!2| , p!3| . 

We study here a prediction of the latter scenarios: each UHECR should point directly to 
its source, which can be at cosmological distances. This was already remarked for the case 
the UHECR primary is a new, neutral, stable or very long lived supersjTumetric hadron of 
mass a few GeV M . It also applies to the mechanism of ref . [O in which neutrinos of energy 



E,y = M^/{2m^) = 4 X 10^"'^(eV/m;^) annihilate with dark matter neutrinos via vv Z 







qq to produce hadronic jets. This is because the opening angle between the incident neutrino 
and a particle of energy E produced in the decay is 59 < Mz/{2E) ~ 10^^. Furthermore 
in both these examples the primary has negligible energy loss except for redshift, so can be 
produced in distant objects such as QSOs. 

By contrast, a proton or nucleus would neither point to an astrophysical source nor be 
associated with a large z QSO, since its scattering from the CBR excessively dissipates its 
energy unless z < 0.01. The rms deflection of a proton of energy E traversing randomly- 
oriented patches of magnetic field having rms value 6B and scale length A is given by 



56 ~ 7.2ydA/200(Mpc)2((55/E)(100EeV/10-^G), (1) 
where d is the distance from the source to the observation point p[. Note that there is some 



evidence for sheets of much higher magnetic field . 

We list in Table | all the UHECR events whose energy is at least 1 a above 8 x 10^^ 
eV and whose direction is known with a solid angle resolution of 10 deg^ or better. The 
energy cut is imposed in order to exclude contamination from events which may be due 
to proton primaries and the angular resolution requirement is necessary to reduce random 
background. The angular resolution in general improves with energy, so both cuts would 
have to be relaxed in order to enlarge the sample. 

A few comments on Table | are in order. The error bars on the energy of the Fly's Eye 
event include systematic as well as statistical uncertainty. The parameters of the Haverah 



Park events are taken from |15|, with errors on the positions determined by us using the 



formulae given in ref. |T^. The angular error in the longitudinal direction, relevant for Afl, 
is Aa*cos6 where Aa is the error in right ascension and 6 is the declination. AGASA reports 
an angular cone radius, denoted cr^ below, defined such that 68% of the events would be 
contained within the error cone: 1° from statistical error alone and 1.6° including systematic 
errors p!7|J^. Our information on AgllO come from ref. 0, which does not give an error on 
the energy measurement, although ref. WE\ quotes a 30% error in general so we expect this 
event satisfies our cut. The high energy event observed by Yakutsk, with E = l.l±0.4x 10^°, 



has too large an angular uncertainty and is not la above 8 x 10 eV so it is not used. See 



| r5| and for other high energy events which we cannot use. 

A correlation with quasars has already been noted for the two highest energy events. 
Elbert and Sommers ^ searched within 10° of the highest energy event, the 320^11 EeV 
event observed by the Fly's Eye group They identified the exceptionally radio- loud 
quasar 3C 147 as an ideal source, aside from its extreme distance. With a redshift of 0.545, 
its distance is of order 2000 Mpc. Biermann |l^ pointed out that another remarkable quasar. 



PG0117+213, is inside the error cone of the second highest energy event (210 EeV) ||17| . 
At 2; = 1.493, about 3500 Mpc, the distance also seemed too great for it to be seriously 
considered as a source. 

The surface density of QSO's is large enough that these two alignments are not statis- 
tically significant and may be accidental. However acceleration of protons to > 10^^ eV 
requires a remarkable source, so if the hypothesis is correct it may be possible to identify 
a more restricted class of sources, with low surface density, for which the correlation is 
statistically significant. 

One of the best-motivated cosmic ray acceleration regions is the jet of an AGN (Active 
Galactic Nucleus), where relativistic shocks and large magnetic fields are found. Depend- 
ing on the age of the AGN, the orientation of its jet with respect to Earth, the "clouds" 
surrounding the inner accretion disk and their relationship to the jets and Earth, and the 
amount of dust in the host galaxy, the same source can be a blazar, radio galaxy, or a quasar. 
It can be unusually bright at visible wavelengths and/or optically variable. The shape of its 
radio spectrum depends on whether it is a full-sized quasar or compact. 

One would like to impose the seemingly trivial criterion that the energy flux in cosmic rays 
implied by the UHECR observation itself, not be much larger than the total electromagnetic 
energy output of the source. However even this simple condition is not straightforward to 
implement for AGNs, due to their directional anisotropy. For instance a blazar pointed away 
from us has a much higher total energy output than evidenced by its observed luminosity. 
Moreover the energy output in a given wavelength band can differ by orders of magnitude 



depending on the intervening material which can "reprocess" the electromagnetic energy. 

In examining the properties of 3C 147 we noticed that it is a compact quasar - that is, 
its jets are only about 1/10 the size of a full-sized quasar with radio lobes. Other indicators 
of its compact character are its optical variability and the fact that its spectrum is cutoff at 



low radio frequencies ||2^ . An anomalous spectrum such as this is characteristic of compact 
radio-loud sources (CSS and GPS) We therefore defined the following specific criteria 
for compact QSOs (CQSOs): 

• QSO in the NASA/IPAC Extragalactic Database (NED). 



• Radio-loud. In practice, we required that the object appear in the Bonn catalog ^1 
This is a compilation containing 1835 radio sources including all those whose flux 
density is > 1 Jy at 5 GHz. The whole sky, excluding the galactic plane < 10°), 
is covered. The surface density of this class of sources is therefore 1835/(34100 deg^) = 
0.054 deg^l 

• Flat or falling radio spectrum at low frequencies. One-third of the Bonn catalog entries 
have a fiat or falling spectrum at low frequencies, so the background surface density 
of the CQSO category is 0.018 deg-^. 

We first determine the probability that the UHECR events actually point directly to 
the candidate sources, given the experimental measurement errors. After that we find the 
probability that randomly distributed compact QSO's, given their surface density, would 
have an equally good alignment to that observed. 

We employ the method of maximum likelihood, which is a standard tool in High Energy 



Physics. For a concise review see the probability and statistics sections of ref. p2|. One 
makes use of the quantity 



where Nu is the total number of UHECR events in the analysis, cr^j is the error on the Xjth 
coordinate, Xi is the measured value of the coordinate (the UHECR position) and x° is the 



(hypothetically) true value of the coordinate, namely the ith source CQSO position. For an 
error cone 0"^ the residual of an event (its contribution to the total x^) is 2.28|r — r°p/cr^. 
The errors on the QSO positions are negligible in comparison with those on the UHECRs. 
A generalization of Eqn. (Q) could be used if correlations in the errors on the coordinates 
of a given UHECR event were available and non- negligible. 

Since there are two degrees of freedom for each event, a residual of about 2 or less corre- 
sponds to a good fit. The expected fluctuations in the sum of the residuals is proportinately 
less than that of any given residual, so that as A*"^ increases the statistical power of the 
analysis increases. For a given set of UHECR events and associated hypothetical sources, 
one determines the Confidence Level of the fit. The Confidence Level (CL) is the probabil- 
ity, with gaussian measurement errors, that an ensemble of = 2Nu measurements will 
produce a large or larger than the observed value. An explicit formula for determining 
the CL corresponding to a given ^ind A"^ is given in the statistics section of p2[. For 
reference, CL = 0.44 for x^ = 10.0 and A'^ = 10, while if x^ = 32.0 and A'^ = 4, CL = 
2 X 10-^ or if = 80.0 and A'^ = 10, CL = 5 x lO^^^. 

Table || gives the residuals {Sx"^) for each of the 5 events listed in Table |I|, under the 
hypothesis the source is the nearest CQSO. As a check of the method, we make the same 
analysis for a second category of "Test" QSOs (TQSOs) chosen to have similar surface 
density and systematics to the CQSOs, by requiring a QSO in NED with 0.400 < z < 0.600. 
This range of z was intentionally chosen to include 3C147, the QSO associated with the Fly's 
Eye event, in order to mimic the CQSO search. By using the same portions of the sky, and 
considering QSO's rather than another type of object, we avoided introducing systematic 
differences between the TQSO and CQSO classes. There are 7 TQSOs in the 5 cones of 
radius 5 deg centered on the 5 UHECR events, giving a surface density of 0.0178 deg~^. Since 
there is no physical motivation that having a redshift in the range 0.4 < z < 0.6 should be 
related to a QSO's acceleration potential, we should NOT find a positive correlation for the 
TQSO category. 



The first row of Table |TTl| gives the probability (CL) to find a total as good as the one 



observed for CQSOs. As a check that the results aren't skewed by having used the properties 
of 3C147 to define the CQSO class, we also give the result when the analysis is restricted to 
the 4 other events. Evidently, the hypothesis that UHECR primaries travel undeflected from 
compact QSOs provides an excellent explanation for the observations, and is equally good 
for the restricted analysis. The same is not true for a randomly chosen category of QSO 
with the same surface density, as evidenced by the very low confidence level (< 2.3 x 10^^) 
for the TQSO fits. The results for this case are shown in the second row of Table |T|. 

By a straightforward Monte Carlo calculation, one can determine the probability distri- 
bution that randomly distributed ohiects having the same surface density as CQSO's, 0.018 
deg~^, produce a given value of x^- The large of the TQSOs is in fact typical of the 
random-background case: the probability to find ^ 61.1 is 0.59. The most interesting 
aspect of the probability distribution is the area below x"^ = 9.02, since this is the prob- 
ability that the CQSO correlation is a statistical fluctuation. The results are given in the 
bottom line of Table fT|. The probability that the correlation observed betwen CQSOs and 
UHECRs is accidental is 0.005Q. Since the correlation hypothesis is a priori, there is no 
reason to restrict to just 4 events; we include that result in Table |T| because the factor-6 
increase in the statistical significance of the observed correlation in going from 4 to 5 events 
gives an idea of the improvement in this analysis due to adding a well-measured UHECR 
event to the sample. 

Let us summarize the underlying assumptions and limitations of the statistical analysis 
presented here. First, we have assumed that the position errors are gaussian and uncorre- 
lated. Therefore our results should be taken as qualitative rather than quantitative indicators 
of the relative probabilities. In the future, cosmic ray experiments should report as detailed 

^Note that the naive procedure of taking the product of the probabilities of finding a random 
source inside each 1-a error region underestimates this probability by several orders of magnitude 
due to neglecting configurations in which some small residuals compensate a large one. 
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information as possible on the positional errors of each high energy event. Secondly, we 
have assumed that the density of compact quasars is approximately uniform. This may not 
be valid, either due to some physical structure or to non-uniformity in the surveys near the 
different UHECRs, although there is no obvious reason to suspect this to be the case. 

It is important to stress that having an incomplete catalog from which to choose the 
best source can only reduce, not exagerate, the quality of the fit if the alignment hypothesis 
is correct. It cannot lead to an incorrect estimate of the random background probability 
as long as the surface density is approximately uniform and is computed from the same 
population as used to find the candidate sources. A dedicated survey within the 3-sigma 
error ellipse of each of the UHECR candidates, and also in several comparable random and 
nearby patches of sky, would be valuable here. 

The candidate sources we have identified should be studied in greater detail and with 
better resolution to learn more about their properties and see if there is a better character- 
ization of the sources. 

The hypothesis that compact radio QSOs are responsible for the highest energy cosmic 
rays gets support from a clustering of events noted by AGASA |0. Three of the five UHECR 
events studied here have one or two companions - nearby events with energy near the GZK 
bound. Near 3C147 we have FE320 at an angular distance of 1.9 deg, Ag62 at 4.0 deg and 
the Yakutsk 120-230 EeV event whose error box is large but overlaps. Near 0109-1-224 we 
have Ag210 at 2.0 deg, HP69 at 4.0 deg and Ag51 at 2.6 deg, and near 1851-1-485 we have 
AgllO at 1.5 deg and Ag43 at 1.5 deg. The lower energy members of these pairs or triplets 
either have a low enough energy to be interpreted as a proton consistent with the GZK 
bound and directions consistent with the angular deflection of protons given in Eqn (P, 
or small enough angular distance from the CQSO source to be interpreted as having been 
undefiected. 

AGASA has announced |]23| the observation of another event with energy above 10^° eV, 



but has not yet released its coordinates or energy. If the correlation pointed out here is 
real, we predict that the new AGASA event will have a compact radio QSO directly behind 
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it, within measurement errors. A radio search in its vicinity may be necessary to check 
this prediction. Since the random probabihty to find a CQSO within a 1 deg cone is 0.05, 
even a few more events with good directional information can confirm or cast doubt on the 
correlation we have found. HiRes will have an angular accuracy 5 times better than Fly's 
Eye and a capacity to collect data which is approximately an order of magnitude higher. 
The Auger detector expects to have a 1 deg resolution, and an even higher data-collection 
rate. 

To summarize, we have found that the highest energy cosmic rays are consistent with 
traveling undefiected from compact radio quasars located at cosmological distances. The 
probabhty that this is a statistical fiuctuation is 0.005. For the moment these results are only 
a tantalizing hint that the highest energy cosmic rays may point directly to their sources and 
travel cosmological distances. However if this hint is born out by future data, nature will have 
revealed some new particle physics mechanism involving neutral, GZK-evading messengers. 
Large statistics with improved angular resolution would allow the properties of the source to 
be more precisely charcterized. UHECRs would then complement traditional astronomical 
tools for studying these extremely distant and powerful sources and their physics. 
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TABLES 



UIIECR 




ENERGY 


RA (dog) 


D(-c (d(-g) 




FE320 


15.10.91 


q 9n+0.92 
"^•^"-'-0.94 


85.2 ±0.5 




2.6 


Ag210 


03.12.93 


(1.7-2.6) 


18.9 


21.1 


8.0 


HP120 


18.04.75 


1.20 ± 0.10 


179± 3 


27± 2.8 


6.7 


AgllO 


06.07.94 


1.10 


280.7 


48.4 


8.0 


HP 105 


12.01.80 


1.05 ± 0.08 


201 ± 8.7 


71 ± 2.5 


7.1 



TABLE 1. Events with E > W'^^eV and solid-angled error less than 10 deg . 



Compact QSO 


Test QSO 


Candidal e 


z 






Candidate^ 


z 


Sep'ii 




3C147 


0.545 


111.6 


1.2 


3C147 


0.545 


111.6 


1.2 


0109-F224 




119.7 


3.5 


0133+207 


0.425 


254.4 


16.3 


1204-F281 


2.177 


138.5 


0.9 


1153+317 


0.418 


286.3 


26.4 


1851+485 


1.25 


89.0 


2.0 


1908+483 


0.513 


254.9 


16.2 


1.345+73 


0.29 


183.4 


1.4 


1300+69 


0.570 


155.1 


0.9 



TABLE II. Compact and Test QSOs nearest the UHECRs of Table I. Separation is in arcmin. 
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Source Class 


xl 


Probability 


xl 


Probability 


Compact QSO 


9.02 


0.53 


7.82 


0.45 


"Test" QSO 


61.1 


2.3 X lO-'^ 


59.9 


4.9 X 10-1° 


Random QSO 


< 9.02 


0.005 


< 7.82 


0.03 



TABLE III. Rows 1,2: Probability for Compact and Test QSOs to produce the observed total 
X^. Row 3: Probability for random sources with surface density of CQSO's (0.018 deg"^) to give 

equal or better than observed. Columns 2,3: all 5 UHECR events; columns 4,5: excluding Fly's 
Eye Event. 
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